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Eighty years ago Eyring proposed that the shear viscosity of a liquid, 77, has a quantum limit rj > nh 
where n is the density of the fluid. Using holographic duality and the AdS/CFT correspondence 
in string theory Kovtun, Son, and Starinets (KSS) conjectured a universal bound ^ for 

the ratio between the shear viscosity and the entropy density, s. Using Dynamical Mean-Field 
Theory (DMFT) we calculate the shear viscosity and entropy density for a fermionic fluid described 
by a single band Hubbard model at half filling. Our calculated shear viscosity as a function of 
temperature is compared with experimental data for liquid ®He. At low temperature the shear 
viscosity is found to be well above the quantum limit and is proportional to the characteristic 
Fermi liquid 1/T^ dependence, where T is the temperature. With increasing temperature and 
interaction strength U there is significant deviation from the Fermi liquid form. Also, the shear 
viscosity violates the quantum limit near the crossover from coherent quasi-particle based transport 
to incoherent transport (the bad metal regime). Finally, the ratio of the shear viscosity to the 
entropy density is found to be comparable to the KSS bound for parameters appropriate to liquid 
®He. However, this bound is found to be strongly violated in the bad metal regime for parameters 
appropriate to lattice electronic systems such as organic charge transfer salts. 

PACS numbers: 71.27.-l-a,05.60.Gg,67.10.Jn 


I. INTRODUCTION 


The viscosity of a fluid is a measure of its resistance 
to externally applied shear or tensile stress. The shear 
viscosity of a fluid measures the resistance of a fluid to 
shear flows, where adjacent layers of a fluid move parallel 
to each other but with different speeds. The differential 
speed between different layers will give rise to friction 
between different layers which will resist their relative 
motion. This is known as the viscous drag. For example, 
the viscous drag force per unit area in the cc-direction, 
Txy, due to velocity gradient dux{y)/dy in the perpen¬ 
dicular j/-direction is given by : 


Txy — rj 


duxjy) 

dy 


( 1 ) 


where rj is the coefficient of shear viscosity. The SI unit 
of shear viscosity is Pascal-seconds (Pa.s) equivalent to 
Newton-second per square meter (N.s m~^). The shear 
viscosity of water is about 10“^(Pa.s) at room tempera¬ 
ture whereas the shear viscosity of highly viscous fluids 
such as glasses near the glass transition temperature can 
be as large at 10^^ Pa.s. 

For fluids rj can be measured through Stokes law for 
sound attenuationji 


a = 


2uj‘^r] 

2,pcl 


( 2 ) 


where a is the rate of attenuation, p is the mass density 
of the fluid, lo and Cg are the frequency and velocity of 
sound in the medium, respectively. This equation has 
been used to determine the shear viscosity as a function of 
temperature for liquid ^He (a correlated neutral fermion 
fluid). Extensive experimental data has been reviewed 
by Huang et ali^. 


The shear viscosity for an electron gas in metals, calcu¬ 
lated from solution of the Boltzmann equation, is given 


rj = -nhkpi, (3) 

5 

where n is the density of electrons, kp is the Fermi veloc¬ 
ity, and f is the electronic mean free path, respectively. In 
the quasi-particle regime of transport kpi ^ 1, i.e., the 
mean free path is much larger than the lattice spacing, 
a ^ kf^. Hence, in analogy with the Mott-Ioffe-Regel 
(MIR) limit, gmir = for minimum metallic con¬ 
ductivity, we can conjecture a lower limit for the shear 
viscosity, pq : 


Pq = ^nh (4) 

corresponding to the case where the electronic mean free 
path becomes comparable to lattice spacing. Also, a com¬ 
parable limit p > nh was proposed by Eyringi almost 80 
years ago. For a large class of strongly correlated sys¬ 
tems like 3d transition metal oxide compounds, organic 
charge transfer salts such as k-(BEDT-TTF) 2X, the MIR 
limit is violated^^— and the coherent quasi-particle based 
transport picture breaks down, i.e. i < a. Similarly, we 
might expect that in the incoherent regime of transport 
the shear viscosity, p, could violate the quantum limit to 
coherent transport, i.e. p < pq. 

Recently a string theory based approach has been 
proposed to understand incoherent quantum transport 
in strongly correlated electron systems, especially the 
strange metal regime of doped cuprates^^—. The key 
idea of this method is to map a strongly coupled confor¬ 
mal field theory (CFT) to weakly coupled gravity in the 
anti-de Sitter (AdS) space in higher dimensioni^. This 
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is known as the holographic duality or AdS/CFT cor¬ 
respondence. Furthermore, event horizon dynamics of 
a black hole in the anti-de Sitter space can be mapped 
to the dynamics of classical fluids. Using the AdS/CFT 
correspondence Kovtun, Son, and Starinets (KSS)^ cal¬ 
culated the ratio, 77 /s, of the shear viscosity ( 77 ) and the 
entropy density (s) in a specific string theory model (type 
IIB) and proposed a universal lower bound for the ratio 


77 h 
- > - 

s dTrfcs 


(5) 


in any material or field theory. This bound is found to be 
well respected in classical fluids like water and quantum 
fluids like the quark-gluon plasma created in the Rela¬ 
tivistic Heavy Ion Collider (RHIC)i^, ultracold atomic 
Fermi gases in the unitary limit of scatteringi^, and by 
theoretical calculations for graphene^^ and for ultracold 
atomic Fermi gase a^^d^ . It has recently been found that 
for the strongly interacting Fermi gas, both experimen- 
tally^o and theoretically 2 ij that the viscosity-entropy ra¬ 
tio is a minimum, not at unitarity, but on the BEC side, 
with a minimum value of ~ 0.2h/kB. Possible violations 
of the KSS bound have been discussed for higher deriva¬ 
tive versions of gravity and with the inclusion of massive 
quarks^^. 

In a recent calculation we tested a jelated but distinct 
bound on charge diffusivity, D > where is the 
Fermi velocity, proposed by Hartnollii. We found^^ clear 
violation of this bound in the strong coupling (bad metal) 
regime of the Hubbard model. In the present paper we 
calculate the shear viscosity in a single band Hubbard 
model and explore possible violations of the conjectured 
quantum bounds on 77 and 77 /s. 

Overall we find that the scale of the viscosity in a cor¬ 
related band system with a lattice constant a in d dimen¬ 
sions is set by 



fmV 

TT 

\mbj 

2d{d-l) 


( 6 ) 


where m is the free fermion mass and mb = h? j{a^Eb) is 
a mass scale determined by some energy scale Eb defined 
by the band structure, such as the half-band width W or 
the rescaled hopping integral t* for a hypercubic lattice. 
A detailed derivation of the expression in Eq.® will be 
provided in the following sections. We show that for a 
lattice system 77 Q is the relevant scale for the analogue of 
the Mott-Ioffe-Regel limit. We will see that the presence 
of this new scale, absent in a conformally invariant sys¬ 
tem, can increase the likelihood of violation of the KSS 
bound. 

The organization of the paper is as follows. In Section 
H we introduce the Kubo formula for calculation of the 
shear viscosity using linear response theory. In Section 
HI we briefly describe the Dynamical Mean-Field Theory 
(DMFT) approach for calculating properties of a single 
band Hubbard model and the iterated perturbation the¬ 
ory (IPT) based approach used to treat the DMFT self- 
consistency for the associated single impurity Anderson 


model. In the same section we introduce calculation of 
the shear viscosity and entropy density in DMFT. In Sec¬ 
tion IV we first briefly review experimental results for the 
temperature dependence of the shear viscosity of liquid 
^He and its possible description by a Hubbard model. In 
Section V we show our results for the Hubbard model on 
the Bethe and hypercubic lattices at half filling. Similar 
results are obtained for both lattices. We compare our 
calculations to experimental results for liquid ^He. The 
temperature dependence of the ratio of the viscosity to 
the entropy density is calculated. It is found that in the 
bad metal regime near the Mott metal-insulator transi¬ 
tion this ratio can be smaller than the KSS bound. In 
Section VI we discuss about experimental measurement 
of shear viscosity in charged systems and finally we con¬ 
clude in Section VH. 


II. SHEAR VISCOSITY 


Nonrelativistic simple fluids are characterized by the 
conserved mass density p, the momentum density tt and 
the energy density £. These quantities will satisfy fol¬ 
lowing conservation laws^i : 


dp 

at 

— V.TT 

(7) 


at 

VjHy 

(8) 

ar. . 

at 

-v.r 

(9) 


where H^ is the momentum current density that the fol¬ 
lowing discussion shows is central to the shear viscosity. 
As a consequence, in analogy with the case of Ohm’s law 
for electrical conductivity : = (JapE^ = —cr„p9p(/(r), 

the generalized Newton’s law for shear flow is 

IIq-p — (fO) 


where rja/s-ys is a viscosity tensor. In particular, the mo¬ 
mentum current density Tlxy in the presence of a trans¬ 
verse velocity gradient is given by^ 


H 


= -It- 


dux 
dy ’ 


( 11 ) 


where 77 = rjxyxy is the coefficient of shear viscosity for 
an isotropic fluid. 

The velocity field Ux{y) gives rise to a perturbation 
with Hamiltonian 

H' = - [ d^r Ux{y)Ttxir) = — [ d^rllxy. (12) 

J too dy J 

To derive Eq. (1121) we have used the conservation law in 
Eq.[51 integration by parts, and 7r(r, t) = exp(—7wt)7r(r). 
The momentum current density IVxy induced by the per¬ 
turbation E[' can be calculated from linear response the¬ 
ory. The shear viscosity is then obtained by taking the 
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limit w —> 0: 

withSi>26 



V 


- lim Im ^ 

uj^O U) 


(13) 


(Pr dt e“‘ 9{t){ (r, t),na;j^(0,0) (14) 


where i/ = is the unit cell volume and 9{t) is the Heav¬ 
iside step function. This formula is the analogue of the 
Kubo expression for the electrical conductivity involving 
the current-current correlation function. 

For a Fermi gas with a quadratic energy dispersion the 
momentum current density operator is given 

^ ( 15 ) 

J m 

where i5/ = /k — /k is the deviation of the distribution 
function from local equilibrium. For Bloch electrons in a 
crystal lattice— in a Bloch state, '0kn(r) = (r|n., k), with 
energy e„(k) 


9en(k) _ 
dka im 


dr tl^lnir) 


dtpjinir) 

dra 


— (n,k|pa|n,k()16) 
m 


where e„(k) is the energy dispersion of the n-th energy 
band. Then (|15l) can be written 


J d'ii Vk^VkJf (17) 

with Vk^ = ^ being the velocity of the Bloch electron. 
Using deformation potential theory^^ a similar result was 
found by Khan and Allen^S when investigating sound at¬ 
tenuation by electrons in metals. 

It should be pointed out that in a general fluid there 
are two terms in the stress energy tensor: one associ¬ 
ated with the kinetic energy and the second with the 
interparticle interaction. In dense classical liquids the 
terms in the Kubo formula due to the interaction term 
dominates and is associated with Einstein-Stokes relation 
where the viscosity is inversely proportional to the parti¬ 
cle self-diffusion constant. In contrast, in dilute gases and 
fluids the kinetic term dominates and the shear viscosity 
scales with the diffusion constant and scattering time^. 
However, for a zero-range interaction, as in the unitary 
Fermi gas (and presumably in the Hubbard model), it 
can be shown that the potential term in the stress tensor 
does not contribute to the shear viscosity.For a discussion 
of the above see around Eq. 7 in Ref. 


III. DYNAMICAL MEAN FIELD THEORY 

We consider the single band Hubbard model with near¬ 
est neighbor hopping, described by the Hamiltonian 

H = -t [cl^Cja + H.C.) - + u'^ ni-fn^lf8) 

{ij),(7 i,o- i 


where Uia = cJ^Cio-, t is the hopping amplitude, p is 
the chemical potential, and U is the Coulomb repulsion 
when a given site is doubly occupied by two fermions 
with opposite spin configuration. Despite its simplicity 
this model has no exact solution except in one dimension. 
The study of this model in higher dimension involves var¬ 
ious approximations. However, as in the case of classical 
mean field theory for the nearest neighbour Ising model, 
in the limit of large dimension, d —>■ oo the model reduces 
to an effective single site model provided we do the scal¬ 
ing t —>■ t* / on a d-dimensional hypercubic lattice^. 

Under this approximation we neglect all spatial fluctua¬ 
tions yet fully retain local quantum dynamics. The self¬ 
energy Ey (u;) for the lattice model then becomes local, 
i.e. Sij(uj) = S(uj)Sij. This is known as the Dynamical 
Mean-Field Theory^ (DMFT) approximation. 

It has been found that DMFT gives a good descrip¬ 
tion of the correlation driven Mott metal-insulator tran¬ 
sition observed in 3d transition metal oxides and the 
crossover from a coherent Fermi liquid to incoherent bad 
metal state with increasing temperature^. Furthermore, 
DMFT has also been found to provide quantitative de¬ 
scription of the resistivity^! and the frequency dependent 
optical conductivity!! for organic charge-transfer salts 
that can be described by a half- filled two-dimensional 
Hubbard model on an anisotropic triangular lattice!!. 
DMFT combined with electronic structure calculations 
based on density functional theory (DFT) has given an 
excellent description of a large class of transition metal 
and rare earth compounds!!. 

The lattice problem under DMFT can be mapped onto 
an effective single impurity Anderson model!! : 

77inip — ^ ~b ^^(k/C^^doo- “t“ H.C.) 

1,(7 1,(7 

-p'^Udoa+ Undotndoi, (19) 


where Udoa = dg^do^. The operators and doo- char¬ 
acterizes a local site and {cl^,cia} characterizes the ef¬ 
fective bath arising from fermions at all other sites. It is 
important to mention that the fictitious bath dispersion 
Cl has no relation to the lattice dispersion, Ck- 

The solution of the impurity problem is the toughest 
part and usually involves use of numerical methods such 
as Quantum Monte Carlo (QMC), exact diagonalization 
(ED), or the numerical renormalization group (NRG). 
We use iterated perturbation theory (IPT)!!!! as it 
is semi-analytical, easy to implement, computationally 
cheap and fast. Yet IPT captures the essential physics in 
the parameter regime U < O.SUc, where Uc is the critical 
value of U at which the zero temperature Mott metal- 
insulator transition happens. Except in close proximity 
of the Mott transition IPT was found to be in good agree¬ 
ment with results from other impurity solvers such as 
the numerical renormalization group (NRG)!! and con¬ 
tinuous time quantum Monte Garlo (CTQMG)!!. In the 
next sub-section we discuss DMFT self- consistency using 
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IPT. 

We briefly mention why it is appropriate to com¬ 
pare the results of AdS/CFT to a calculation Involv¬ 
ing DMFT. The latter becomes exact in infinite dimen¬ 
sion. Generally, AdS/CFT Is concerned with finite- 
dimensional quantum field theories. However, It Is found 
that in certain parameter regimes DMFT can accurately 
give a quantitative description of quasi-two-dimensional 
metals near the Mott insulator. Furthermore, connec¬ 
tions have been made between the results of AdS/CFT 
and the infinite-dimensional limit of a model for a gap¬ 
less spin liquid^. We know that the DMFT approxima¬ 
tion reduces a lattice problem to an effective local Im¬ 
purity problem which captures local correlation effects. 
On the other hand the AdS/CFT correspondence maps 
a strongly coupled field theory to a problem of fluid me¬ 
chanics and fluids are characterised by short range cor¬ 
relations. So we might expect DMFT based description 
of quantum transport of lattice electronic systems will 
be closely related to quantum transport in the hydrody¬ 
namic regime of a strongly coupled held theory. 


symmetric case, which we consider In the present study, 
(iv) The new self-energy can be calculated using IPT 
ansata^ as 


S(a;) = Un + 


AE( 2 )(a;) 

1 - HE( 2 )(u;) 


(23) 


where 

^ n(l-n) ^ D(l-n)-/r + /ro 

no(l-no)’ no(l-no)C/2 ^ ^ 

and n, no are the local and bath particle numbers, re¬ 
spectively. Is the self energy from second order 

perturbation theory and Is given by 


+00 2 

f l[{de,po{e,)) 


nF(-ei)nF(e2)n_F(-e3) 
uj + i 6 — ei + 62 — €3 


^ nF(ei)nF(-e2)nF(e3) 

a; -I- — ei -I- £2 — £3 


(25) 


A. Iterated perturbation theory 


The irreducible self-energy In IPT Is approximated us¬ 
ing the second order (in U) polarization bubble involving 
fully interacting bath Green’s function Go(w). The self¬ 
energy under this approximation can be shown (using 
moment expansion of the interacting density of states) 
to smoothly interpolate between the atomic limit t = 0 
and the weak-coupling limit t/ —?► 0. In the following 
paragraph we briehy discuss DMFT self-consistency us¬ 
ing IPT as the impurity solver. As we are Interested in 
calculating transport properties we work with real fre¬ 
quencies, as against the imaginary frequency formula¬ 
tion that requires analytical continuation of Imaginary 
frequency data to real frequency. 

(1) For a given lattice density of states iVo(£) and self¬ 
energy E(u;) the local Green’s function is given by 


G(w) 


^+00 


No{e)de 


a;+ + — e — S(a;+)' 


( 20 ) 


where fi is the local chemical potential and = uj + iS 
with (5 > 0. 

(ii) From the knowledge of the local Green’s function, 
G(u;), we can calculate the bath hybridization function, 
A(w) by using 

A(a;) = + fJ- — S(w) — G~^(a;) (21) 


(iii) Subsequently using bath hybridization we can calcu¬ 
late bath Green’s function as 


Go(w) 


I 

UJ + po — A(a;) 


( 22 ) 


The parameter /xg = p, — Un is the bath chemical po¬ 
tential and It vanishes at half-filling for the particle-hole 


where po{uj) = — ilm[Go(a;+)] and (5 —^ 0+. We iterate 
(i)-(iv) until the desired self-consistency in self-energy 
and other physical quantities are achieved. Here we 
consider the particle-hole symmetric case at half filling 
n = 1. In this case p = ^ ioi all U and T. 


B. Shear viscosity in DMFT 


Using the self-consistent self-energy we can calculate 
the shear viscosity. In the limit of d —> 00 all vertex cor¬ 
rections to two-body correlation functions drops oul4^ 
and the temperature dependent coefficient of shear vis¬ 
cosity, rjlT), given by the Kubo formula Eq. (|T^ can be 
calculated using a simple polarization bubble as 


rj{T) = — r°°duj 

^ J —00 


dnp(uj) 

duj 


f*+oo 


de Qa;y{e)A'^{uj,el^ 6 ) 


where f = is the volume of the unit cell of a d- 
dimensional hypercubic lattice with lattice constant a. 


A{w, e) = -Im 


I 


_uj+ + p - e - S(a;+)_ 


npiuj) = 


1 


(27) 

(28) 


g/ 3 c^ -p I 

are the spectral density and Fermi function, respectively. 


Qxy{e) = “ Ek) 


(29) 


with Vka = is the transport density of states for 

the shear viscosity and N is the number of lattice sites. 
Following a similar procedure to that in Ref. we can 
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show that the transport density of states for shear viscos¬ 
ity for a d- dimensional hypercubic lattice with nearest 
neighbour hopping is given by 


Qxyif) — 


7 


2 d{d-l) 


--M3(e) + 2eM2(e) 


-I- (4c?t^ —-e^)Mi(e) — 4t^eMo(e) 


where 7 = and 


(30) 


e 

M„(e)= J z^No{z)dz (31) 


where No{e) = ~ ^k) is the density of states per 

spin. In the Appendix we give a detailed derivation of 
this important result. In the following sub-sections we 
explicitly evaluate this expression for the hypercubic lat¬ 
tice and Bethe lattice cases. 

One should consider how the vertex corrections could 
modify the DMFT results in finite dimensions. For the 
unitary fermi gas vertex corrections increase the viscosity 
by a factor of about 2.6 [compare the discussion below 
equation (54) in Reference |3l|. For the quark-gluon plas- 
mon in the theory of Quantum Chromodynamics (QCD) 
at high temperatures vertex corrections significantly in¬ 
crease the viscosity, changing the functional dependence 
on the coupling constant [compare equations (4.25) and 
(4.26) in Reference . In a two-dimensional Fermi liq¬ 
uid the vertex corrections have been shown^ to be of the 
order of {\n[EFIT)Y. For a doped Hubbard model it was 
found in a Dynamical Cluster Approximation calculation 
based on a four site cluster that the vertex corrections to 
the optical conductivity were not significant, except very 
close to the Mott insulator—. A study of the same model 
using a two-particle self-consistent approach found that 
vertex corrections changed the calculated resistivity by 
less than a factor of two^. In light of the above it seems 
unlikely that vertex corrections would increase the vis¬ 
cosity by more than an order of magnitude compared to 
the DMFT results. 


1. Hypercubic lattice case 


are the dimensionless integrals and i = ejt* is the dimen¬ 
sionless energy. 

We define the scaled dimensionless transport density 
of states for the shear viscosity, 0 (e), as : 


0 a:i/(e) — 

Using the 


-3/3(?) + 4e'J2(?) - - 2)/i(g) - 2?/o(e')(34) 


exact integrals 



[ z^e-^"dz = -l-{e^ + l)e-^" 
J — 00 ^ 


(35) 

(36) 

(37) 


we get 


0 xy(Q 




(38) 


In Fig. [T] (a) we show transport density of states for vis¬ 
cosity for hypercubic lattice. Interestingly, the transport 
density of states for electrical conductivity for hypercu¬ 
bic lattice also follows a relation similar to Eq. (1551) . as 
shown in Ref. 15 

The shear viscosity is then given by 


^+00 


r] = Vo 


doj 


dnpidj) 

dCj 


r+oo 


de Qxyie)A‘^{uj,eX39) 


where the dimension full prefactor tjq is given by Eq. ([5]) 
with nib = and ui = uj/t* is the dimensionless energy. 


2. Bethe lattice case 

We consider the Bethe lattice (Cayley tree) with coor¬ 
dination number z. In the limit of infinite coordination 
number (z 00 ), the density of states has semicircular 
form^ : 


Eor hypercubic lattice in the limit of d —>■ 00 we have 
A^o(e) = exp[—for the density of states. It is 
important to mention that the chosen density of states 
in the limit d —>■ 00 requires the scaling : 2 t t*/\fd. 
The transport density of states for the shear viscosity is 
then given by 

Qxy {() = [-373(e) + 47/2(7) - 7^/1 (7) 

+2/1 (e) — 27/0(7)] (32) 


where 


In{x) = 



(33) 


No{^) = ^Vw^e{W-\e\), (40) 

where 6 {x) is the Heaviside step function, W = 2t* is 
the half bandwidth and the nearest neighbour hopping 
amplitude (t) in this case is scaled as t —>■ t* Eor a 

Bethe lattice with coordination number z the connectiv¬ 
ity K = z — 1 while that for a d-dimensional hypercubic 
lattice is 2d. So, in the limit of large coordination num¬ 
ber we can always do the mapping z O 2d. Because of 
its tree like structure the Bethe lattice has no closed loop 
and hence no energy dispersion with Bloch wavevector k. 
However, by invoking the /-sum rule we can still calcu¬ 
late 0 x 1 /(e). For the given density of states we then have 


M^e"" du 
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the following exact integral 


,44 


= ie7Vo(e) + i +-tan ^ 

A Z TT 

Mi{e) = -^{W^ - e^)Noie) 

M 2 (e) = —^---^7Vo(e) + — 


_VW 2 -e\ 


(41) 

(42) 



Maie) 


—— tan 
47r 


e 


1 


No(e) [l8e'‘ - 5e^W^ 


16iy^] 


(43) 

(44) 


Then by replacing these exact analytical integrals into 
the expression in Eq. (1501) for Oxyie) and using W = 
2 t\/^ we get 



Oxyis) 


24iM(j3T)«o(e)[8ir‘-25.“»'^ + 26.*](45) 


It is interesting to mention that the constant term as 


well as the tan 




term cancels out in the final 


expression for Oxyie). 

In Fig. [I] (b) we show the scaled dimensionless trans¬ 
port density of states, Oxyie) '■ 


Oxy(e) 


1 

IW 


Noie) {8 - 25e^ + 26i^) 


(46) 


FIG. 1: (Color online) Density of states, No{e), and scaled 
dimensionless transport density of states for shear viscosity, 
Qxy{e), for the hypercubic lattice [panel (a)] and for the Bethe 
lattice [panel (b)]. Oxyie) for the Bethe lattice shows addi¬ 
tional structures near the band edges but for the hypercubic 
lattice it is just proportional to the density of states. Oxyie) 
for the Bethe lattice has been multiplied by a factor of 10 
to show both curves in the same panels. Energies have been 
scaled by the effective hopping amplitude, t*, in the case of 
the hypercubic lattice and by the half-band width, W, in the 
case of the Bethe lattice. 


with e = tjW for the Bethe lattice. Near the band edges 
(e = ±1) Oxyie) shows non-monotonic structures in con¬ 
trast to the density of states, A^o(e) = f a/I — e^, which 
is always monotonic near the band edges. 

The expression in Eq. (1261) for the shear viscosity for 
the Bethe lattice is then given by Eq. (l39l) where the 
dimensionfull prefactor tjq is given by Eq. (jS]) with rrii, = 
and w = uj/W is the dimensionless energy. 


C. Entropy density 


The total internal energy in DMFT is given by^ 


EjT) 

N 


+ 00 


kBTj2 / de 


eNoie) 


iuj„ + ^^(iWn) - e 


+ -'^T.aii0Jn)Gcriiu}n) (47) 


where N is the total number of particles in the system, 
LOn = (2n + Vj'KksT is the Matsubara frequency, and 
A^o(e) is the non-interacting density of states. In the 
paramagnetic state Eq. dUl) can be expressed as a real 
frequency integral 


EjT) 

N 


/ +ca 

dujnpiu!)iLjj + ^)Aioj) 

-oo 


TT 


(iwnF(w)Im [E(a;~'’)G(a;~'')] , (48) 


where Aiuj) = — :^Ini [G(w“'')] is the spectral function. 
From E(T) we can calculate the specific heat using 

CviT) = ^ ^ ^ and then we can calculate the local 

entropy density, s(T), as 


siT) 


1 

V 



CxiT') 

T' 


dT', 


(49) 


where v is the volume of the system. The temperature 
dependence of the specific heat and the entropy for the 
half-filled Hubbard model have both been calculated pre¬ 
viously using a range of impurity solvers including IPT 
and Quantum Monte Carlo^^— . It is found that in the 
metallic phase the entropy density is linear in tempera¬ 
ture below the Fermi liquid coherence temperature, Tcoh 
and becomes of order nks ln(2) for T ~ Tcoh, where n is 
number density of fermions. 


D. Quantum limits 

The quantum limit of the shear viscosity, = ^nh, 

1 2 

is based on the free particle dispersion Ek = in the 
continuum limit. For a discrete lattice model we need to 
derive an appropriate quantum limit for shear viscosity. 

For temperatures and frequencies much less than the 
coherence scale (i.e. T <C Tcoh, w ksTcoh where Tcoh 
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is the coherence temperature which is of the order of the 
Kondo temperature for the corresponding single impurity 
Anderson model) the self energy, S(a;), has the Fermi 
liquid form : 

I](w, T)=uj(^l--iC [uj^ + {TTkeT)^] , (50) 

where Z is the quasi-particle renormalization factor and 
C is a positive constant. 

Following the procedure in Ref. d used to estimate the 
Mott-Ioffe-Regel limit for the conductivity we can show 
that at low temperature (T <C Tcoh) the shear viscosity 
for the hypercubic lattice is given by 

r,{T) = r]^oQxy{0)Ioi7Th*T{T)/h (51) 

where /oi ~ 0 .08 is a dimensionless integral and 

h/riT) = —Im(i;(w = 0,T)) is the quasi-particle de¬ 
cay rate. The quantum limit to shear viscosity will then 
correspond to t* ^ h/T{T) and we will have the quantum 
limit to shear viscosity 

for the hypercubic lattice and 

< = (53) 

for the Bethe lattice. 


IV. PARAMETERS FOR LIQUID ®He 

We consider liquid ^He because of the availability of ex¬ 
tensive experimental data for the temperature and pres¬ 
sure dependence of the shear viscosity, recently reviewed 
and parametrised by Huang et ah— . First, we review how 
liquid ^He might be described as a lattice gas with a 
Hubbard model Hamiltonian. 

Low temperature properties of liquid ^He can be de¬ 
scribed by Landau‘s Fermi liquid theory. The effective 
mass of the quasi-particles [as deduced from the spe¬ 
cific heat] is about 3 times the bare mass m at 0 bar 
pressure and increases to 6 times at 33 bar, when the 
liquid becomes solid. The compressibility is also renor¬ 
malised and decreases significantly with increasing pres¬ 
sure. This led Anderson and Brinkman to propose that 
^He was an “almost localised” Fermi liquid. Thirty years 
ago, Vollhardt worked this idea out in detail, consider¬ 
ing how these properties might be described by a lattice 
gas model with a Hubbard Hamiltoniani^ The system 
is at half filling with U increasing with pressure, and 
the solidification transition (complete localisation of the 
fermions) then has some connection to the Mott transi¬ 
tion. All of the calculations of Vollhardt were at the level 
of the Gutzwiller approximation (equivalent to Kotliar- 
Ruckenstein slave boson mean-field theory). A signifi¬ 
cant result from the theory is that it describes the weak 


pressure dependence and value of the Sommerfeld-Wilson 
ratio of the spin susceptibility to the specific heat [which 
is related to the Fermi liquid parameter Fq]. At ambient 
pressure U was estimated to about 80 per cent of the crit¬ 
ical value Uc for the Mott transition. Vollhardt, Wolfle, 
and Anderson^ also considered a more realistic situation 
where the system is not at half-filling. Then, the doping 
(band filling) is determined by the ratio of the molar 
volume of the liquid to the molar volume of the solid 
(which by definition corresponds to half filling). Finite 
temperatures extension to Volhardt theory was done by 
Seiler, Gros, Rice, Ueda, and Vollhardti^ Later Georges 
and Laloux^ argued ^He is a Mott-Stoner liquid, i.e., 
one also needs to take into account the exchange interac¬ 
tion and proximity to a Stoner ferromagnetic instability. 
If this Mott-Hubbard picture is valid for ^He then one 
should also see a crossover from a Fermi liquid to a “bad 
metal” with increasing temperature. Specifically, above 
some “coherence” temperature Tcoh, the quasi-particle 
picture breaks down. For example, the specific heat per 
atom should increase linearly with temperature up to a 
value of order around Tcoh, and then decrease with in¬ 
creasing temperature. Indeed one does see this crossover 
in experimental data (compare Figure 1 in Ref. 

We now consider what Hubbard model parameters 
are appropriate for ^He. The density at a pressure of 
1 bar, n ~ Na/{57 cm^) (where Na = 6.023 x 10^^ 
is the Avagadros number) increases monotonically to 
n ~ Na/{‘26 cm^) at 33 bar (near the solidification pres¬ 
sure) (see Table HI in Ref. [5^ . 

The band mass nib can be written in terms of Ep = 
h^kp/(2m), the non-interacting Fermi energy, and the 
band energy Eb as 


m 

mb 


1/ 7 \2 

e; 




(54) 


where we have used the fact that n = 1/a^ = k^/iiir^). 

There are several ways to estimate the band energy 
scale. If we have a Bethe lattice, then Eb = W = Ep, 
at half filling. Alternatively, we can compare the non¬ 
interacting density of states per spin at the Fermi en¬ 
ergy No{0). This has the value of 3/(4 Af), 2 /(7rIF), 
and l/{y/Trt*) for the cases of a parabolic band (free 
fermions), Bethe lattice, and hypercubic lattice, respec¬ 
tively. Setting these equal gives Eb = W = SEp/^Sn) ~ 
0.85Ep and Eb = t* = AEp/ {3^/tt) ~ 0.75Ep. Using 
the density at I bar and the non-interacting expression 
Ep = h^k^/{2m) we estimate Tp ~ 4.95 K, and so 
t* ~ 3.72 K. 

In the following sections we compare some of our calcu¬ 
lations of the shear viscosity with experimental data for 
^He. Huang et ah— showed that the shear viscosity of sat¬ 
urated liquid ^He from 3 mK to 0.1 K follows the Fermi 
liquid relation rj ex IjT'^. Furthermore, they showed that 
the shear viscosity data in the range from 3 mK to near 
the critical point at 3.31 K, collected over the past 50 
years from various experimental groups can be fitted to 



the empirical form : 


= |4 + ^ + ^+C4 (55) 

with Cl = 2.897 x 10“^ Pa.sK^, C 2 = —7.02 x 10“^ 
Pa.sK^ ®, C 3 = 2.012 x 10"® Pa.sK and C 4 = 1.323 x 10"® 
Pa.s. We note that at low temperatures Eq. (1551) has 
a Fermi liquid term. At high temperatures Eq. (1551) 
has the asymptotic value of C 4 which is comparable to 
nfi ~ 1.6 X 10“® Pa.s at 1 bar pressure. It should be 
pointed out that this is for data along the liquid-vapour 
curve and so the pressure gradually increases with tem¬ 
perature. However, as the critical pressure is about 100 
kPa, much less than the melting pressure, this pressure 
dependence is not significant. The viscosity decreases by 
a factor of at most ten as the pressure increases from 
1 kPa to 3 MPa [the melting pressure] for all tempera¬ 
tures below 1 K. Huang et al. fitted all the avaialable 
experimental data to an expansion in terms of Cheby- 
shev polynomials and used this to plot the temperature 
dependence for pressures ranging from 1 kPa to 20 MPa 
(compare Figure 8 in Ref. Q)- For pressures larger than 
about 500 kPa, the viscosity has a non-monotonic tem¬ 
perature dependence with a minimum around a temper¬ 
ature of 1 K. 


V. RESULTS 

We consider the case of half-filling, n = 1, i.e. each 
site on the average is occupied by one fermion. We study 
the shear viscosity and the entropy density as a function 
of correlation strength, 17, and temperature, T. 

In Fig. [2] (a) and Fig. |2Kb) we show the scaled shear 
viscosity rjiT )/tjq as a function of temperature for various 
interaction strengths U, for the hypercubic and the Bethe 
lattice, respectively. Similar results are obtained for both 
lattices. 


A. Quantum limits 

We consider violation of quantum limit of shear vis¬ 
cosity. In the weakly correlated hypercubic lattice sys¬ 
tem with U = 0.5, the shear viscosity is always above 
quantum limit, ~ 0 . 28 r 7 § but as we increase the in¬ 
teraction strength, U, the shear viscosity smoothly goes 
below the quantum limit with increasing temperature, T. 
This corresponds to the fact that at low temperatures 
(T <g; Tcoh) the quantum transport is due to coherent 
quasi-particle states but at high temperatures {T > Tcoh) 
the transport becomes incoherent in nature. This is the 
analogue of how in bad metals the resistivity smoothly 
increases above the Mott-Ioffe-Regel limit. 




FIG. 2: (Color online) Shear viscosity, ri{T)/rjQ, as a function 
of temperature for a range of interaction strengths, U. Results 
are shown for both the hypercubic lattice [panel (a)] and the 
Bethe lattice [panel (b)[. Solid lines corresponds to DMFT 
based numerical results. T and U are measured in units of 
W for the Bethe lattice and in units of t* for the hypercubic 
lattice case. Note that as the Mott insulator is approached 
the shear viscosity becomes extremely small. The dashed lines 
are the quantum limits given by Eq. (1521) [hypercubic lattice] 
and Eq. (1531) [Bethe lattice]. 


B. Low temperature behaviour 

Fig. [3] clearly shows that the shear viscosity follows 
Fermi liquid characteristic 1/T^ behaviour in the low 
temperature region (T <C Tcoh)- The range of Fermi 
liquid behaviour decreases with increasing U. This is be¬ 
cause the coherence scale [and Kondo temperature for 
the corresponding single impurity Anderson model] de¬ 
creases with increasing correlation strength U. The 1/T^ 
behaviour is similar to the low temperature behaviour of 
the electrical conductivity and the quantum transport 
in this region can be characterised by coherent quasi¬ 
particle states. 

Our calculated shear viscosity shows qualitative be¬ 
haviour consistent with experimental data for liquid 
®He, using parameters estimated in the previous section. 
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There is qualitative as well as quantitative agreement at 
low temperatures but not at higher temperatures. In¬ 
terestingly, our calculated shear viscosity for [/ = 2 for 
the hypercubic lattice nearly fits with the experimental 
results at low temperatures. Our calculation suggests 
^He is a moderately correlated system with U/Uc ~ 0.5 
\Uc ^ 4.0 for the hypercubic lattice] as against the sug¬ 
gestion of Volhardli^ that ^He is a nearly localized Fermi 
liquid with U/Uc ~ 0.8 at 1 bar pressure and U/Uc 0.9 
close to the melting pressure. It is important to men¬ 
tion that Gutzwiller based static mean field theory over 
estimates local correlation effects but the self-consistent 
treatment of dynamic correlation effects in DMFT renor¬ 
malizes local correlation effects. 


C. High temperature behaviour 

In the high temperature region, T ^ Tcoh, the shear 
viscosity shows significant deviation from the low tem¬ 
perature Fermi liquid behaviour as can be observed from 
Fig. [2| and Fig. [3| The quantum transport in this region 
is incoherent in nature. For the weakly and moderately 
correlated systems the deviation is smooth and mono¬ 
tonic but for strongly correlated systems for U = 2.5 and 
above the deviation is much sharper and non-monotonic. 
This is due to the sharp crossover between the Fermi 
liquid fixed point and the local moment fixed point in 
the strongly correlated regime. A similar non-monotonic 
temperature dependence is seen in the electrical resistiv¬ 
ity from DMFT calculations and in organic charge charge 
transfer salts close to the Mott insulator . 



D. Entropy density 

In Fig. m we show the entropy density, s(T), as a func¬ 
tion of temperature for various interaction strengths. At 
high temperatures the entropy density approaches ln(4) 
which arises due to local charge and spin fluctuations. As 
the temperature decreases charge fluctuations freeze out 
and the model can be described by localised weakly inter¬ 
acting spin I/2s with characteristic entropy density ln(2). 
Finally in the Fermi liquid state the local spin degrees of 
freedom are dynamically screened and the entropy den¬ 
sity vanishes linearly in temperature. For weakly and 
moderately correlated electron system the entropy den¬ 
sity smoothly crosses over ln(2). But for strongly corre¬ 
lated electron systems with U = 2.5 and above a kink 
like feature develops. This corresponds to formation of 
poorly screened local moment. The position of the kink 
in the specific heat versus temperature curve is related 
to the coherence temperature, Tcoh^ For extremely cor¬ 
related systems with U = 3.0 and above the entropy 
density given by iterated perturbation theory (IPT) is 
under estimated. Consequently the specific heat in the 
coherent-incoherent crossover region becomes negative, 
which is unphysical. This is due to an incorrect total 
energy estimate in IPT which has been reported in ear¬ 
lier literature. (See for example. Figure 7 in Ref. 

In the unphysical temperature range we set the specific 
heat to zero and the calculated entropy density which 
is an integrated quantity will deviate by not more than 
5% from the actual value. Such a small error has little 
effect on whether the Kovtun-Son-Starinet (KSS) bound 
is violated. 



FIG. 3: (Color online) Shear viscosity, rj{T), on a log-log plot 
clearly shows 1/T^ behaviour in the low temperature region, 
characteristic of a Fermi liquid. At high temperatures there 
is significant deviation from the Fermi liquid behaviour as 
the Mott transition is approached. The upper dashed line 
corresponds to experimental data for saturated liquid ^He 
parametrised by Eq. ([55]). The lower dashed line is at a 
higher pressure. The results are shown for the hypercubic lat¬ 
tice case with the parameters t* = 3.72 K and m(,/m = 3.6. 
Both T and U are measured in units of t*. 


FIG. 4: (Color online) Entropy density, s{T), (in units of 
nks) as a function of temperature, T, for various interaction 
strength, U. Below the coherence temperature, Tcoh, the en¬ 
tropy is linear in temperature, characteristic of a Fermi liquid. 
The crucial point is that for T > T^ok, the entropy is of order 
nks- The kink like feature for U = 2.5 and above corresponds 
to the formation of poorly screened local moment and its po¬ 
sition is closely related to Tcoh- The calculation is for the 
hypercubic lattice and both T and U are measured in units 
of t*. 
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E. Possible violation of the KSS bound 


Finally we consider the dimensionless scaled shear vis¬ 
cosity, r]{T), entropy density, s(T), ratio 




(56) 


At the KSS bound ({T) = (for d = 3). As stressed 

before this depends on the material properties t* and a 
as well as the mass of the fermion, m. 

In Fig. Owe show ({T) and compare to its value against 
the Kovtun-Son-Starinet (KSS) limit for parameters ap¬ 
propriate for ^He and typical lattice electronic systems 
like cuprates and organic superconductors. 

For cuprates^ the hopping integral t ~ O.lSeV, a = 
3.9A and for organic charge transfer salts^ t ~ O.OSeV, 
a = SA. For these systems d = 2 and hence t* = 2t^fd = 
2\/2t. This will give mjrrH) ~ 1.0 for cuprates and 
m/rrib — 1.2 for organics, as compared to m/m}, ~ 3.6 
for liquid ^He. As a result the shear viscosity for these 
lattice electronic systems will be smaller by a factor of 
about 10 than for the charge neutral fermionic fluid ^He. 
Straub and Harrison considered a simple model for the 
hopping integral for d-bands in transition metal^ This 
gives for d-sigma bands m/mb — ‘2-S{rd/a)^, where 
is approximately the d-state radius and of the order of 
the radius transition metal atom, ~ lA. In principle 
then for a system with a large lattice constant, the band 
width can be very small and values of m/mb even smaller 
than unity are possible. In an ultracold fermionic atom 
system in an optical lattice one could in principle then 
make m/mb, and thus the viscosity, extremely small. 

From Fig. [5] we can clearly see that for all 17 <3.0 
and for ^He parameters C,{T) is above the KSS limit. For 
extremely correlated system U = 3.5 there is strong vio¬ 
lation of the limit in the crossover region but even for this 
system at high temperature the bound seems to be re¬ 
spected (within numerical error in calculation of entropy 
density, s{T)). Also, in the high temperature region the 
scaled ratio seems to approach some universal limit. 

For electronic lattice systems the limit is well respected 
in the coherent quasi-particle regime of transport but the 
limit is violated in the region T > T^oh- This is due to re¬ 
duction of the shear viscosity by a factor of 10 compared 
to ^He parameters. The violation is as large as 1000% 
for these systems, when they are close to the Mott tran¬ 
sitions. 


VI. EXPERIMENTAL DETERMINATION OF r? 

IN ELECTRONIC SYSTEMS 

Given our result that the KSS bound can be violated 
in a bad metal it is highly desirable that experimental 
measurements be performed on candidate strongly cor¬ 
related electron materials such as organic charge trans¬ 
fer salts and cuprates. Unfortunately, at present there 



FIG. 5: (Color online) Dimensionless scaled ratio, C,(T), of 
the shear viscosity and entropy density. [See Eq. (I56II for 
definition of C(7')]. It is shown as a function of temperature, 
T, for various interaction strength, U. Dashed and dot-dashed 
lines corresponds to the Kovtun-Son-Starinets (KSS) limit for 
^He and typical lattice electronic systems, respectively. For 
electronic systems there is strong violation of the quantum 
bound for T Tcoh- The calculation is for the hypercubic 
lattice case and both T and U are measured in units of t*. 


is no direct measurement of the shear viscosity for elec¬ 
tronic systems. Recently an indirect estimate of 77 /s was 
made from angle resolved photo-emission spectroscopy 
(ARPES) experiments in cuprates^ giving a value com¬ 
parable to the KSS limit. However, it should be stressed 
that neither the viscosity nor the entropy were directly 
measured. Rather, the ARPES lineshape was used to es¬ 
timate the quasi-particle lifetime and the state occupa¬ 
tion. The viscosity was then estimated from the lifetime. 
The entropy was estimated from an expression in terms 
of the state occupations in a non-interacting fermion sys¬ 
tem. In the incoherent regime of transport this method 
will not be applicable. It is important to mention that 
our calculation showed that in the coherent quasi-particle 
regime of transport 77 /s is always above the KSS bound 
and hence our result is consistent with these experimen¬ 
tal results. 

A more direct way to measure the viscosity of the elec¬ 
tron fluid in a metallic crystal is through the attenuation 
of sound, as first emphasized by Mason^. A more sophis¬ 
ticated and general theory was developed by Kahn and 
Alienas. The connection between shear viscosity and ul¬ 
trasound attenuation can be loosely motivated by Stokes 
law, given in Eq. Q. In a metal, provided the wave¬ 
length of sound is much larger than the electronic mean 
free path, then one is in the hydrodynamic limit, and 
the attenuation is given by a similar expression to Stokes 
law (with appropriate indices for crystal axes), with p the 
solid density, not that of the electron fluid. In a simple 
free electron model Eq. ([3]) shows that that the electronic 
viscosity is proportional to the scattering time, just like 
the electrical conductivity. Hence, the ultrasound at- 
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tenuation should scale with the conductivity. Indirect 
evidence for this idea was found from the temperature 
dependence of ultrasound attenuation in aluminiunt^, 
including the predicted quadratic frequency dependence. 
In clean metals the attenuation (and viscosity) becomes 
very large at low temperatures, making it easier to mea¬ 
sure. Also, for high frequency ultrasound, one can reach 
the “quantum regime” where the mean free path becomes 
comparable to the sound wavelength. Pippard worked 
out a general theory describing the crossover from the 
hydrodynamic regime to this quantum regime^. In bad 
metals could one experimentally measure the small elec¬ 
tronic viscosity, of the order of nhl First, the small mean 
free path, characteristic of bad metals, means one will al¬ 
ways be in the hydrodynamic regime. However, the small 
viscosity means that the sound attenuation due to the 
electron fluid will be small and possibly dominated by 
other sources of attenuation such as crystal dislocations. 
A rough estimate for an electron viscosity of order of nh, 
and a sound frequency of 1 GHz gives an attenuation of 
less than 0.1 cm“^, of the order of typical sensitivity, such 
as in measurements for heavy fermion compounds^. 

Resonant Ultrasound Spectroscopy (RUS)^“— has 
been used to make measurements on strongly correlated 
electron systems^^— . The spectrum is determined by 
the resonant elastic modes of the sample; they are deter¬ 
mined by the sample shape and orientation, elastic con¬ 
stants, and dissipation. RUS allows determination of the 
elastic constant tensor Cij from measurements on small 
samples (< 1 mm^ volume). In the regime where the 
attenuation of the ultrasound is dominated by coupling 
to the electrons, rather than fluctuations associated with 
phase transitions, the viscosity could be determined from 
the damping [frequency width, Aw] of the resonances. 
We estimate Auj/uj ^ rjuj/Cij^ and so for the MHz fre¬ 
quencies typically used in RUS, the damping associated 
with a viscosity of order nh requires an oscillator Q factor 
of order 10^° and so is unlikely to be observable. 

Recently, several new approaches have been suggested 
to experimentally measure the viscosity of the electron 
fluid in a metallic crystal. Forcella, Zaanen, Valentinis, 
and van Der Marel^i considered electromagnetic proper¬ 
ties of viscous charged fluids, finding signatures due to 
the viscosity such as negative refraction, a frequency de¬ 
pendent peak in the reflection coefficient, and a strong 
frequency dependence of the phase. However, they note 
that these effects may be difficult to observe for viscosi¬ 
ties of the order of nh. Tomadin, Vignale, and Polintf^ 
considered a two-dimensional electron fluid in a Corbino 
disk device in the presence of an oscillating magnetic flux. 
They showed that the viscosity could be determined from 
the dc potential difference that arises between the in¬ 
ner and the outer edge of the disk. In particular, for 
viscosities of the order of nh the potential difference var¬ 
ied significantly oscillation frequencies in the MHz range. 
Levitov and Falkovicb^S recently considered the flow of 
an electron fluid in a micrometer scale channel in the hy¬ 
drodynamic regime, where the electron-electron collision 


rate is much larger than the momentum relaxation rate. 
They found that when the viscosity to resistance ratio 
is sufficiently large viscous flow occurs producing vortic- 
ity and a negative nonlocal voltage. Spatially resolved 
measurements of the voltage allow determination of the 
viscosity. Torre, Tomadin, Geim, and Polini22 consid¬ 
ered the electron liquid in graphene in the hydrodynamic 
regime and showed that the shear viscosity could be de¬ 
termined from measurements of non-local resistances in 
multi-terminal Hall bar devices. Although these propos¬ 
als are promising for the two-dimensional electron fluids 
in graphene and semiconductor heterostructures fabrica¬ 
tion of the relevant micron-scale devices may be partic¬ 
ularly challenging for bad metals such as cuprates and 
organic charge transfer salts. 

VII. CONCLUSIONS 

We have studied the shear viscosity, the entropy den¬ 
sity and their ratio for a single band Hubbard model 
using single site Dynamical Mean Field Theory. Sim¬ 
ilar results were obtained for the density of states as¬ 
sociated with both hypercubic and Bethe lattices. We 
compared our results for the temperature dependence of 
the shear viscosity to experimental results for liquid ^He. 
The calculated shear viscosity shows qualitative as well as 
quantitative behaviour consistent with experimental re¬ 
sults. At low temperatures the shear viscosity is propor¬ 
tional to 1/T^ corresponding to coherent quasi-particle 
based transport in the Fermi liquid state. At high tem¬ 
peratures the shear viscosity shows significant deviation 
from Fermi liquid state behavior. This corresponds to 
crossover from coherent quasi-particle based transport 
to incoherent transport (the “bad metal”). With in¬ 
creasing interaction strength U the shear viscosity be¬ 
comes less than conjectured quantum limits of shear vis¬ 
cosity, of the order of nh. Finally, we considered the 
scaled dimensionless ratio between shear viscosity and 
entropy density. This ratio in the Hubbard model de¬ 
pends on the energy scale t*, length scale a, and the 
free fermion mass m. This is in contrast to the universal 
limit predicted by Kovtun, Son, and Starinets using 
the AdS/CFT correspondence in a conformally symmet¬ 
ric field theory model. For ^He parameters the ratio is 
above the universal bound but for parameters appropri¬ 
ate for electronic lattice systems, such as cuprate and or¬ 
ganic metals, this bound is found to be strongly violated, 
in the bad metal regime near the Mott metal-insulator 
transition. We hope that our results will stimulate ex¬ 
perimental measurements of the shear viscosity in bad 
metals. 
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Using the convolution theorem we can easily show that 
each term of Yiuj) has the following form 

+ 00 

= J Fa{e - z)Ga{z)dz, a = I,--- ,3(63) 

— OO 

where. 


APPENDIX: Transport density of states for shear 
viscosity 

In the limit d —> oo the viscosity will involve the fol¬ 
lowing transport function 

&xy{e) = m^'^vlvlS{e-ek). (57) 

k 

For a d-dimensional hypercubic lattice 

2 /Oi '|4 4 

^xy{e) = - ^4 - ^ sin^ikx) sm'^{ky)S{e - Ck). (58) 

k 

To evaluate this we first we define the Fourier transform 


E„(e) = T F„(cc)e*‘"^dw 


' —OO 

p+oo 


Ga(e) = / G„(a;)e“^da;, 


(64) 


and Ya{Lu) = Fa{ui)Ga{ui). 
For the hrst term 


+00 

= ^ / i 




= sgn{e) =-^sgn{e) (65) 


+ 00 


r(w) = J 0,y(e)e-“^de 

— OO 

d 

= 72(2t)4^sin2(fei)sin2(fc2) Y[ 

k 

Jl{2tLu) 


= 7^(2t)^Jg ^{ 2 tuj) 

(2^7)^ 4.-2 


CX — 1 
1 2 


2 tuj 


J^-\2tuj) [Ji{2tuj)Y 


(59) 


where 7 = ^ 


Using relations for Bessel functions 

1 d?J^{2tu)) 


+00 

If dVoi 2 tu)]\^^ 

^ ^ 


— OO 
.2 1 


= -e^Noie) 
where we have used 

d[Jo(2ta;)]^ 


( 66 ) 


duj 


+00 


= 0 


[Jo(2H]''e-^|:“=0. 


(67) 


Finally we get 


[Jo{2tu;)f duj^ &Y^ie) = J z^No{z)dz 


-f—— —jQ{ 2 toj) — 


d- 1 °2t(d- l)a; 
we can rewrite U(w) as 


^ ^ [Jo(2M]^”^ Ji(2t4i^) 


(2t7)^e 

4(d-l) 


d(d — 1) 


z'^No{z)dz, (68) 


~ d{d- 1) ^ + Id - Fr z^No{z)dz = J2i^el = 2t^d and 


^ [M2tuj)\^-^ Ji{2tu:) 


(d-1): 


= Yi{uj)+Y2{uj)+Y3{uj). (61) 

We can Fourier transform back to calculate 0xy(e) as 

+00 

0.y(e) = ^ / y{^)e^‘^"duj. (62) 


/_r z^No{z)dz = 0. 

A similar exercise for the second term will give 

eg'(.) = 


(d-1) 

(2^7)^ . 

(d-l) J 


2(d- 1) 

O 

zNo{z)dz. (69) 
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For the third term we have 


Finally, we have 


+ 00 


= lim 77 - • 7 ri • sgn{e) 

uj^o 2 -k 2! dw2 L \ ^ ^ > 


(70) 


and 


G3(e) = - 


1 2^7^ 

27r [d — 1) 

1 7' 


+ 00 


[Jo(2ta;)]‘^ ^Ji(2tw)da; 


+00 

^ f d[Jo(2ta;)]'^ 


27r d{d — 1) 


duj 


duj 


7 


d{d-l) 


eiVo(e). 


(71) 



f 

2d{d-l) J 


zNo{z)dz 


(2t7)^£ 

4(d-l) 


+ 




d{d-l) 


e 

J z'^No{z)dz — 


2d{d- 1) 


z^No{z)di72) 


Collecting and rearranging all the terms we obtain Eq. 

dsni). 
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